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Abstract: The behavior of hadronic matter at high baryon densities is studied within 
^ 1 Ultrarelativistic Quantum Molecular Dynamics (URQMD). Baryonic stopping is observed 

for Au+Au collisions from SIS up to SPS energies. The excitation function of flow shows 
O . strong sensitivities to the underlying equation of state (EOS), allowing for systematic 

studies of the EOS. Dilepton spectra are calculated with and without shifting the p pole. 
^ [ Except for S+Au collisions our calculations reproduce the CERES data. 
G\ '• 

p3 1 Introduction 

The only possibility to probe excited nuclear matter in the laboratory are nucleus-nucleus 
reactions [Q. In particular when two heavy ions like Au or Pb collide most centrally, the 
combined system forms a zone of high (energy) density and high excitation of the involved 
constituents. The transient pressure at high density has specific dynamic implications, 
such as collective sideward flow. Hence, fundamental properties like the repulsion of the 
nuclear equation of state (EOS) are studied via event shape analysis of nucleons and 
clusters The EOS at fixed temperature yields a density dependent potential and a 

modified nucleon mass. At low densities these effects are proposed by chiral lagrangians||. 
|7|]. Since the chiral condensate (qq) relates closely to hadron masses, the decay of short 
lived vector mesons, observed through the dilepton channel, is suggested as a promising 
experimental signal to investigate the gradual restoration of chiral symmetry. 

2 Ultrarelativistic Quantum Molecular Dynamics 

Since many important aspects of nuclear matter are not observable, numerical transport 
models are suited to test which assumptions are compatible to nature. The present model 
(URQMD) includes explicitely 50 different baryon species (nucleon, delta, hyperon 
and their resonances up to masses of 2.11 GeV) and 25 different meson species (including 
strange meson resonances), which are supplemented by all isosp in-projected states (see 
Table 1). Symmetries regarding time inversion, iso-spin, charge conjugation, etc. are 
implemented in a general manner, e.g. all corresponding antiparticles are included and 
treated on the very same (charge-conjugate) footing. For excitations of higher masses 
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Table 1: List of implemented baryons, mesons and 
their resonances. In addition all charge conjugate 
and iso-spin projected states (and photons) are 
taken in and treated on the same footing. 
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Figure 1: 02 partial decay rates into 
specific channels. The average life- 
time is given by the inverse of the 
sum. Hence, in URQMD particles 
below resonance mass live longer, 
due to shrinking phase space. 




a newly developed string model is invoked. It consistently allows for the population 
of all included hadrons from a decaying string. At low energies the dominant part of 
MM and MB interactions is modeled via s-channel reactions (formation and decays of 
resonances), whereas BB interactions are designed as exchange of charge, strangeness 
and four momentum in the t-channel. For all resonances we use mass-dependent decay 
widths as illustrated in Fig.|I] for the 02 meson. The lifetime of resonances is calculated 
as their inverse width. There exist, however, recent theoretical ansatzes which yield a 
different mass dependence for the life-times of resonances |K|. The real part of the 
baryon optical potential is modeled according to the simple Skyrme ansatz, including 
Yukawa and Coulomb forces. 

3 Creation of dense nuclear matter: stopping 

Baryonic stopping is a necessary condition for the creation of hot and dense nuclear matter. 
The key observable is the rapidity distribution of baryons. It is displayed in Fig.|2] and 
[| for heavy systems such as Au+Au and Pb+Pb at energies referring to three presently 
used heavy ion accelerators. In all cases gaussian rapidity distributions with peak around 
midrapidity are found. However, the physical processes associated show characteristic 
differences: The average longitudinal momentum loss in the SIS energy regime is mainly 
due to the creation of transverse momentum, whereas at AGS/SPS energies abundant 
particle production consumes a considerable amount of the incident beam energy. 

At CERN/SPS energies baryon stopping is influenced also by the formation time of 
strings which are excited in hard collisions. In URQMD baryons originating from a leading 
constituent (di-)quark at the string edges interact with (2/3)1/3 and mesons with 1/2 of 
their full cross sections during their formation time r. The sensitivity on this reduction 
is shown in Fig.|3| for the system S+S at 200 AGeV. The default calculation (including 
formation time) reproduces the data |TT[ fairly well whereas the calculation with zero 
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Figure 2: Rapidity distributions for 
Au+Au collisions at SIS (1 AGeV), AGS 
(10.6 AGeV) and Pb+Pb at CERN/SPS 
energies (160 AGeV). 
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Figure 3: Rapidity distribution for 
S(200 AGeV)S for various treatments of 
the constituent (di-)quark cross section 
(see text). 



formation time (dotted line) exhibits strongest stopping. A calculation with zero cross 
section within the formation time gives transparency. 

In order to study the influence of this this effect more closely the \fs distributions for 
Au+Au collisions at AGS and S+S collisions at SPS energies are analyzed. Fig.|5| (right) 
shows the respective distribution for Au+Au. The collision spectrum is dominated by BB 
collisons with full cross sections and exhibits a maximum at low energies. Approximately 
20% of the collisions involve a diquark, i.e. a baryon originating from a string decay whose 
cross section is reduced to 2/3 of its full cross section. 
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Figure 4: Rapidity distributions for S+S at 200 AGeV (left) and Pb+Pb at 160 AGeV 
(right). The histograms label from top to bottom: negative hadrons {h~), kaons (K~ , 
K + ), protons (b + — b~) and lambdas (A). The kaons are multiplied by five. The symbols 



show data from NA35 

In Fig.|5] (left) the same analysis is performed for S+S at 200 AGeV. In contrast to the 
heavy system at AGS the collision spectrum exhibits two pronounced peaks dominated 
by full BB collisions, one in the beam energy range and one in the low (thermal) energy 
range. Now approximately 50% of the collisions, most of them at intermediate i/i values, 
involve baryons stemming from string excitations whose cross sections are reduced by 
factors of 2/3 (referred to as diquarks) or 1/3 (referred to as quarks). The peak at high 
y/s values stems from the initial hard collisions whereas the peak at low energies is related 
to the late, thermal stages of the reaction. 
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Figure 5: ^fs distributions for baryon baryon collisions in central reactions of Au+Au 
(left) and S+S (right) at AGS and SPS energies respectively. 

4 Probing the repulsion of the EOS: flow 

The creation of transverse flow is strongly correlated to the underlying EOS [|TJ. In par- 
ticular it is believed that secondary minima as well as the quark-hadron phase transition 
lead to a weakening of the collective sideward flow. The occurrence of a phase transition 
should therefore be observable through abnormal behaviour (e.g. jumps) of the strength 
of collective motion of the matter fi~3f . Note that URQMD in its present form does not 
include any phase transition explicitly. In Fig.^| the averaged in plane transverse momen- 
tum is displayed for Au+Au from 0.1 to 4 AGeV incident kinetic energy. Calculations 
employing a hard EOS (full squares) are compared to cascade simulations (full circles). 
In the latter case only a slight energy dependence is observed. In contrast, the calcula- 
tion with a hard EOS shows strong sensitivity. Here, the integrated directed transverse 
momentum per nucleon is more than twice as high as for the cascade calculation. This 
indicates the importance of a non-trivial equation of state of hadronic matter. 

The amount of directed transverse momentum scales in the same way as the total 
transverse momentum produced in the course of the reaction. Hence, the directivity 
depends only on the reaction geometry but not on the incident energy. This is demon- 
strateted in Fig. [7|, where the mean p x as a function of the rapidity divided by the average 
transverse momentum of all particles is plotted. 
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Figure 6: Excitation function of the total 
directed transverse momentum transfer px- 
dir for Au+Au. URQMD calculations in- 
cluding a hard EOS (full squares) are com- 
pared to the predictions of cascade calcu- 
lations (full circles). 
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Figure 7: Mean directed transverse mo- 
mentum as a function of the scaled rapid- 
ity. The transverse flow < p x > (y/y pro j.) 
scales with the mean transverse momentum 
< p t >, i.e. directivity does not depend on 
the bombarding energy. 



5 Temperature dependence of the EOS: photons 

Semiclassical cascade models in terms of scattering hadrons have proven to be rather 
accurate in explaining experimental data. Therefore it is of fundamental interest to extract 
the equation of state from such a microscopic model, i.e. to investigate the equilibrium 
limits and bulk properties, which are not an explicit input to the non-equilibrium transport 
approach with its complicated collision term (unlike e.g. in hydrodynamics |T3|,[I3). In 
Fig||the thermodynamic properties of infinite nuclear matter are studied within URQMD. 

Infinite hadronic matter is simulated in URQMD by constructing a box of 250 fm 3 vol- 
ume with periodic boundary conditions. According to the saturation density, nucleons are 
initialized randomly in phase space, such that a given energy density is reproduced. After 
the system has equilibrated according to the simulation with URQMD the temperature 
is extracted by fitting the particles' momentum spectra. Alternatively, the temperature 
can be extracted from the relative abundances of different hadrons, e.g. the A/N ratio. 

In Fig.|8] the result of this procedure is compared to various analytic forms of the EOS. 
While the EOS of a Hagedorn gas and a QGP yields energy densities e ~ lGeV/fm 3 at 
T = 150 MeV the temperature dependence is much smaller in URQMD. It yields about 4- 
5 times less energy density, being in fair agreement with a gas composed of nonrelativistic 
nucleons and ultrarelativistic pions. It remains to be seen whether a reparametrization of 
the resonance continuum in the Hagedorn model as suggested in Ref. [|I7[] would resolve the 
deviation as compared to URQMD. On the other hand, beyond T ~ 200 MeV the energy 
density rises much faster than T 4 approaching even the QGP value of e ~ 10 GeV/fm 3 
around T = 300 MeV. This indicates an increase in the number of degrees of freedom. 
It may be interpreted as a consequence of the numerous high mass resonances and string 
excitations, which seem to release constituent quark degrees of freedom (but, of course, no 
free current quarks as in an ideal QGP). Investigations of equilibration times and relative 
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Figure 8: 'EOS' of infinite nuclear matter 
function of the energy density versus 
temperature at fixed net-baryon density of 
p B = 0.16/fm 3 in URQMD (symbols). The 
curves refer to analytical forms of the EOS, 
i.e. a Hagedorn-gas (top), a quark-gluon 
plasma (middle), and an ideal gas of nu- 
cleons and ultrarelativistic pions (bottom). 
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Figure 9: Transverse momentum spectrum 
of directly produced photons in Pb+Pb 
collisions at 160 AGeV calculated with 
URQMD. The resulting spectrum is com- 
pared with hydrodynamical calculations. 
In all models the processes ni] \— > rcy, 
up i— > 7T7 and tttt i— > p7 are considered as 
photon sources. 



particle and cluster abundances are in progress. Moreover, the admittedly poor statistics 
have to be improved, in order to study the high temperature behavior. 



Experimentally, the EOS can be accessed by measuring electromagnetic radiation JT5 
In Fig.y the direct photon production from meson+meson collisions in Pb+Pb collisions at 
160 AGeV is shown. Here, only mesons stemming from string decays are included. Elastic 
meson-meson scattering with a d = 15mb (independent of y/s) was allowed. The result 
is compared to calculations within the 3-fluid model jlljl , scaling and Landau expansion 
with Ti = 300 MeV. 

6 In medium masses: dileptons 

In Fig.|TU] and [TT] calculations of dilepton spectra with URQMD are shown for p+Be and 
S+Au. Dilepton sources considered here are Dalitz decays (jr , t] and u) and vector 
meson decays (p, uo and 0). Dalitz decays of heavier meson and baryon resonances are 
included explicitely via their emission of p mesons (assuming vector meson dominance). 
In order to avoid double counting, the p mesons from 77's, and uu's are excluded from the p 
contribution. Pion annihilation is included dynamically into the contribution of decaying 
p mesons (tt + tt~ i— > p 1— > e + e~). 

While the result for p+Be agrees well with the published data from CERES/SPS ||18|| , 
two points around M ~ 400 MeV are missed by about two standard deviations for S+Au. 
Speculations about the origin of this deviation include electromagnetic bremsstrahlung, 
annihilations of pions and a modification of the p meson propagator due to a gradual 
restoration of the chiral symmetry. 
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Figure 10: Dilepton mass spectrum for 
p+Be at 450 GeV. The calculation in- 
cludes Dalitz decays and conversion of vec- 
tor mesons (see also legend for S+Au). The 
sum of all contributions (solid curve) is 
folded with the CERES mass resolution. 
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Figure 11: Dilepton mass spectrum for 
S+Au at 200 AGeV (see also legend for 
p+Be). Here no in- medium modifications 
of the p propagator is considered. Around 
M ~ 400 MeV two points are missed by 
< 2a. 
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In URQMD the contribution of pion annihilation to the p-peak 
40% for S+Au. Major additional sources are decays of heavy baryons (A*/N* 
proposed in Ref.|19| and meson resonances (see also Fig.|I]): 
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In Ref.0,0] a linear dependence of the p°/uj pole position as a function of the nuclear 
density p has been suggested: m p o(g/Qo) = m p o(0)(l — Xq/qq). Here p denotes the 
ground state density of nuclear matter, and A = 0.18, in agreement with various other 
calculations. Since the restriction to low densities may not be suitable for heavy ion 
collisions, the following extrapolation towards higher densities is taken: 



m p o(g/g ) 



m p o (0) 
+ A q/qq 



(2) 



In Fig.|12| an application of Eq.@ is made to calculate a dielectron mass spectrum for a 
density dependent vector meson pole. This result yields only a small enhancement around 
M ~ 500 MeV as compared to the calculation without pole shift (bottom curve). On the 
other hand, the data can nicely be reproduced, if the strong (unphysical) assumption is 
made, that the pole at the decay point (p \— > e + e~) is shifted according to the creation 
density (upper curve). This would be a neglection of the finite decay length. The dis- 
crepancy of a calculation without decay length (Ad ec = fm) as compared to the result 
including the decay length is driven by two reasons: i) The increase of the p lifetime 
(~ 7 fm/c) below its resonance mass in the region M ~ 0.3 — 0.5 GeV (where a dilepton 
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Figure 12: Dilepton mass spectrum for 
S+Au at 200 AGeV. The curves label sim- 
ulations with pole shift according to the 
creation density (top), the decay density 
(middle) and without pole shift (bottom). 
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Figure 13: Dependency of p-interactions 
on the density expectation values. The 
curves show the values for the creation, ab- 
sorption and decay of p mesons. Note that 
(g b ) < 0.5g around M ~ 0.4 GeV. 



excess in S+Au is reported |]T8|) lowers the decay density down to (g) ~ 0.2g o for S+Au 
(or 0.3g o for Pb+Au). ii) An enhancement of the decay length leads to an increase of 
reabsorption. Hence, the radiation path for p i— > e + e~ is substantially truncated. This 
fact is further investigated in FigJfJ], where the effect of the mass dependent p-width 
is depicted by the course of the density dependence for p mesons with invariant masses 
M. The decreasing width (i.e. increasing lifetime) of low-mass-resonances leads to higher 
mean decay times where the baryon density is already dilute. Thus in-medium-corrections 
can only yield small enhancements when treated this way. However, the interpretation of 
the experimental data gives the following impression: The data for light systems such as 



p+Be and p+Au (see also Ref . pllJ22]] ) as well as the data for the heavy Pb+Au system 
both for inclusive and central reactions are reproducible without mass shifts. In contrast, 
the central data for S+Au exceed the URQMD calculation around M ~ 0.4 GeV by about 
two standard deviations. 

In URQMD the p meson pole position is shifted according to the density at which the 
p meson decays, i.e. eventually converts into e + e~ . Note that this procedure is equivalent 
to a "shining" description, where the p constantly emits e + e~ pairs according to the rate 
dN ee /dt = r(p i — > e + e~). In Fig.H the shining description is compared to the treatment 
in URQMD for a p source at rest including a time dependent reabsorbation probability. 
Both methods yield - whithin the statistical limits - the same average emission times. 
Furthermore a simple time dependent density profile is used to calculate the mean e + e~ 
emission density. Again, both methods yield the same values. 

The results for Pb+Au are shown in Figs.|I^ and IB]. Both calculations for inclusive 
reactions are in fair agreement with the preliminary observation from CERES |2(J. The re- 
sult for central events is given in Fig.|TS|. Note that the calculations without modifications 
are well compatible for p+Be and Pb+Au for both centralities. Only in S+Au reactions 
two datapoints are missed by the default calculation by about two standard deviations. 
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Figure 14: Comparision of the 
treatment of shining p-sources to 
the URQMD approach, where the 
dileptons are created at the p de- 
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Both methods show the same time 
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7 Summary 

Studies of the equation of state and consequences of gradual restoration of the chiral sym- 
metry are presented using a novel microscopic phase space model, URQMD, including 75 
hadron species and strings. The directed tranverse momentum shows strong sensitivities 
to the underlying EOS: It is small in the cascade calculation, whereas it scales linearly 
with the average transverse momentum for a hard equation of state. Hence, measureing 
the excitation function of the transverse directed flow allows for a systematic study of the 
EOS. The calculation of dilepton yields without modifications of the p mass pole is well 
compatible with the CERES-data for p+Be and Pb+Au. Only in S+Au reactions two 
datapoints at 0.3 — 0.5 GeV are missed by the default URQMD-calculation by about two 
standard deviations. 
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